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ABSTRACT 

We present follow-up spectroscopy and photometry of 11 post common enve- 
lope binary (PCEB) candidates identified from multiple Sloan Digital Sky Survey 
(SDSS) spectroscopy in an earlier paper. Radial velocity measurements using the 
Nai AA 8183.27,8194.81 absorption doublet were performed for nine of these systems 
and provided measurements of six orbital periods in the range P or b = 2.7 — 17.4 h. 
Three PCEB candidates did not show significant radial velocity variations in the 
follow-up data, and we discuss the implications for the use of SDSS spectroscopy 
alone to identify PCEBs. Differential photometry confirmed one of our spectroscopic 
orbital periods and provided one additional P or b measurement. Binary parameters 
are estimated for the seven objects for which we have measured the orbital period 
and the radial velocity amplitude of the low-mass companion star, A" scc . So far, we 
have published nine SDSS PCEBs orbital periods, all of them P or b < 1 d. We perform 
Monte-Carlo simulations and show that 3ct SDSS radial velocity variations should 
still be detectable for systems in the orbital period range of P or b ~ 1 — 10 days. 
Consequently, our results suggest that the number of PCEBs decreases considerably 
for P or b > 1 day, and that during the common envelope phase the orbital energy of 
the binary star is maybe less efficiently used to expell the envelope than frequently 
assumed. 

Key words: Binaries: spectroscopic - stars:low-mass - stars: white dwarfs - binaries: 
close - stars: post-AGB - stars: evolution variables 



1 INTRODUCTION 

Most of the stellar sources are formed as parts of binary or 
multiple systems. Therefore the study of binary star evolu- 
tion represents an important part of studying stellar evolu- 
tion. While the majority of the wide main sequence binaries 
evolve as if they were single stars and never interact, a small 
fraction are believed to undergo mass transfer interactions. 
Once the more massive main sequence star becomes a red 
giant it eventually overfills its Roche lobe. Dynamically un- 
stable mass transfer exceeding the Eddington limit ensues 
onto the companion star, which consequently also overfills 
its own Roche lobe. T he two stars then orbit inside a com- 
mon envelope (CE, e.g. |PaczvnskJll97fj : |l7vio fc Sokerlll988l: 



llben fc Liviol Il993l : iTaam fc Rickerl 120061 ; IWebbinkl 120071 ). 
and friction inside this envelope causes a rapid decrease 
of the binary separation. Orbital energy and angular mo- 
mentum are extracted from the binary orbit and lead to 
the ejection of the envelope, exposing a post common enve- 
lope binary (PCEB). After the envelope is expelled PCEBs 
keep on evolving towards shorter orbital periods through 
angular momentum loss via gravitational radiation and, for 
companion stars above the fully convective mass limit, by 
magnetic wind braking. In binaries that do not undergo a 
CE phase both components evolve almost like single stars 
and keep their wide binary separations. A predicted conse- 
quence is hence a strongly bi-modal binary separation and 
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orbital period distribution among post-main sequence bina- 
ries, with PCEBs concentrated at short orbital periods, and 
non-PCEBs at long orbital periods. 

The scenario outlined above is thought to be a funda- 
mental formation channel for a wide range of astronomical 
objects such as low-mass X-ray binaries, double degener- 
ate white dwarfs, neutron star binaries, cataclysmic vari- 
ables and super-soft X-ray sources. Some of these objects 
will eventually end their lives as type la supernovae and 
short gamma-ray bursts, which are of great importance for 
cosmological studies. 

While the basic scenario of CE evolutio n has been 
known for a long time (e.g. IPaczvnskH Il976h . the phys- 
ical details involved are complex and still poorly under- 
stood. As a result, full hydrodynamical models for the 
CE phase h ave been calculated only for a small num - 
ber of cases |Sandouist et al l l200d : iRicker fc Taamll2008l) . 
Similarly, orbital evolution through magnetic braking is 
not well understood, with different prescriptions in an- 
gular rnranentuin_Joss differ i ng by orders of mag nitudes 
(e.g. IVerbunt fc Zwaanl Il98ll ; lAndronov et al.l [2003V Con- 
sequently, current binary population synthesis models still 
have to rely on simple parametrisations of energy or angular 
mom e ntum equations |Nelemans et al.ll2000l; Dewi fc TaurisI 



2007) 



200d; iNelemans fc Tout! 120051 : iPolitano fc Weilerl |20o3 



A fundamental problem in advancing our understanding 
of CE evolution and magnetic braking in close binaries is the 
shortage of stringent observational constraints that can be 
used to test and calibrate the theory. PCEBs consisting of a 
white dwarf and a main sequence star are a very well-suited 
class of objects to provide the required innovative observa- 
tional input since they are numerous, well understood in 
terms of their stellar components, they are bright and hence 
accessible with 2-8 metre tele scopes, and their study is not 
complicated by mass transfer. ISchreiber fc Gansickd (|2003l ) 
analysed a sample of 30 well studied PCEBs and showed that 
the population of PCEBs known so far is not only small but 
also heavily biased towards young systems with low-mass 
secondary stars, and can hen ce not be used for com parison 
with population models (e.g. IWillems fc Kolb|[2004T ). 
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The Sloan Digita l Sky Survey (lYork et ail 
IStoughton et ail |2002| ; lAdelman-McCarthv et ail 
SDSS) is currently dramatically increasing the number of 
known white dwarf plus main sequence ( WDMS) binaries 
dSilvestri et al 1 l2007l ; ISchreiber et al.|[2007l) , paving the way 
for large-scale observational PCEB population studies. We 
have initiated an observational programme to identify the 
PCEBs among the SDSS WDMS binaries, and to determine 
their binary parame ters (|Rebassa-Mansergas et ail 120071; 
ISchreiber et ail 120081 ) . Those systems in which significant 
radial velocity variation is detected, are classified as close bi- 
naries or PCEBs. Here we present follow-up observations of 
11 PC EB candidates identified in iRebassa-Mansergas et ail 
(2007) and provide accurate values of the orbital periods 
as well as estimates of their stellar parameters and orbital 
inclinations for seven of these systems. 



2 OBSERVATIONS 

We have obtained time-resolved spectroscopy and pho- 
tometry for 11 SDSS WDMS binaries (TableQ] Fig.©, 
henceforth designated SDSS J0052-0053, SDSS J0246+0041, 
SDSSJ0309-0101, SDSSJ0820+4314, SDSS J0314-0111, 
SDSSJ1138-0011, SDSSJ1151-0007, SDSS J1529+0020, 
SDSSJ1724+5620, SDSS J2241+0027 and SDSS J2339- 
0020. Below we briefly describe the instrumentation used 
for the observations, and outline the reduction of the data. 



2.1 Spectroscopy 

(i) Magellan-Baade. Intermediate-resolution long-slit spec- 
troscopy of SDSSJ0314-0111 and SDSS J2241+0027 was 
obtained on the nights of 2007 October 2 and 3 using 
the IMACS imaging spectrograph attached to the Magel- 
lan Baade telescope at Las Campanas Observatory. The 
600^ mm -1 red-sensitive grating was used along with the 
slit- view camera and a slit width of 0.75 arcsec, giving a re- 
ciprocal dispersion of 0.39 A px -1 . The IMACS detector is 
a mosaic of eight 2kx4k SITe CCDs, and long-slit spectra 
with this instrument are spread over the short axis of four of 
these CCDs. Using a grating tilt of 14.7° allowed us to posi- 
tion the NaiAA 8183.27,8194.81 doublet towards the centre 
of CCD2 and the Ha line near the centre of CCD4. 

The images were reduced using STARLINK so ftware, and 
the PAMELA and molly packages l|Marshl [l989h were em- 
ployed to optimally extract and calibrate the spectra. By fit- 
ting Gaussian functions to arc and sky emission lines we find 
that the spectra have a resolution of approximately 1.6 A. 
Helium arc lamp exposures were taken at the start of each 
night in order to derive a wavel ength solution with a statisti- 
cal uncertainty of only 0.002 A (|Marsh et aill 19941 ). This was 
applied to each spectrum taken on the same night. Detector 
flexure was measured and removed from the wavelength so- 
lution fo£_eacli_s_rje£tnrm usjri^ the positions of sky em ission 
lines (ISouthworth et ail 12009 : ISchreiber et ail |200S| ). Flux 
calibration and telluric line removal was performed using 
molly and spectra of the standard star BD +28° 4211. 

(ii) Magellan- Clay and Very Large Telescope (VLT). 
The observations of SDSS J1138-0011, SDSS J1151-0007 
and SDSS J1529+0020 were obtained with the same 
setup and reduced in the same manner as described by 
ISchreiber et ail (|2008r ). and we refer the reader to that paper 
for full details. 

(iii) New Technology Telescope (NTT). Two observing 
runs were carried out at the NTT in August and Octo- 
ber 2007, providing intermediate time-resolved spectroscopy 
of SDSS J0052-0053, SDSS J0246+0041, SDSS J0309-0101, 
SDSS J2241+0027 and SDSS J2339-0020. We used the 
EMMI spectrograph equipped with the Grat#7 grating, the 
MIT/LL red mosaic detector, and a 1 arcsec wide long-slit, 
resulting in a wavelength coverage of A7770 - 8830 A. The 
data were bias-corrected and flat-fielded using the STARLINK 
software, and spectra were optimally extracted using the 
PAMELA package. Helium- Agon arc lamp spectra were taken 
at the beginning of each night and were then used to estab- 
lish a generic pixel-wavelength relation using molly. Specif- 
ically, we fitted a fourth order polynomial which gave rms 
smaller than 0.02A for all spectra. We then used the night 
sky emission lines to adjust the zero-point of the wavelength 
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Table 1. Log of the observations. Included are the target names, the SDSS ugriz psf-magnitudes, the period of the observations, the 
telescope/instrument setup, the exposure time, and the number of exposures. 
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Notes: M-Baade and M-Clay refer to the two Magellan telescopes at Las Campanas observatory. We use sm to indicate that the data 
were taken in service mode. CA2.2 is the 2.2 metre telescope at Calar Alto observatory. 



calibration, thus correcting for instrument flexure. The spec- 
tral resolution of our instrumental setup determined from 
the sky lines is 2.8A. Finally, the spectra were calibrated 
and corrected for telluric absorption within molly using ob- 
servations of the standard star Feige 110. 

(iv) William Herschel Telescope (WHT). One spectrum 
of SDSS J1138-0011 was taken with the 4.2 m William Her- 
schel Telescope in June 2007 at the Roque de los Muchachos 
observatory on La Palma. The double-beam ISIS spectro- 
graph was equipped with the R158R and the R300B grat- 
ings, and a larcsec long-slit, providing a wavelength cov- 
erage of A7600 — 9000 A. The spectral resolution measured 
from the sky lines is 1.6 A. Reduction and calibration were 
carried in the same way as described for the NTT above. 

(v) SDSS. For one system, SDSS J1724+5620, we 
were able to determine an accurate orbital period from 
our photometry alone, but were lacking follow-up spec- 
troscopy. SDSS DR6 contains three 1-d calibrated spec- 
tra for SDSS J1724+5620 (MJD-PLT-FIB 51813-357-579, 
51818-358-318, and 51997-367-564). Each of these spectra 
is combined from at least 3 individual exposures of 900 s, 
which are also individually released as part of DRrJ]. For 
SDSS J1724+5620, a total of 23 sub-spectra are available, of 
which 22 allowed reliable radial velocity measurements. 



berg (Germany) for a total of 13 nights during May, August 
and September 2006 and March 2007. The IAC80 cm tele- 
scope was equipped with a 2k x 2k CCD. A binning factor of 
2 was applied in both spatial directions and only a region 
of 270x270 (binned) pixels of size 0.66 arcsec was read. The 
detector used at the Babelsberg 70 cm telescope was a cryo- 
genic lkxlk Tek CCD. The whole frame was read with a 
binning factor of 3, resulting in a scale of 1.41 arcsec/pixel. 
A semiautomated pipeline involving DoPHOT was used to 
reduce the images and extract the photometric information. 

(ii) Calar Alto 2.2m telescope. We used CAFOS with 
the SITe 2kx2k pixel CCD camera on the 2.2 metre tele- 
scope at the Calar Alto observatory to obtain filter-less dif- 
ferential photometry of SDSS J0314-0111. Only a small part 
of the CCD was read out in order to improve the time reso- 
lu tion. The data were re duced using the pipeline described 
in iGansicke et all ([2004) , which pre-processes the raw im- 
ages in MIDA S and extracts aperture p hotometry using the 
Sextractor (|Bertin fc Arnoutslll996f ). 

(iii) Kryoneri 1.2 m telescope Filter-less photometry of 
SDSS J0820+4311 was obtained in November 2006 at the 
1.2 m Kryoneri telescope using a Photometries SI-502 516 x 
516 pixel camera. The data reduction was carried out in the 
same way as described above for the Calar Alto observations. 



2.2 Photometry 

(i) IAC80 and AIP 70 cm telescopes. We obtained differential 
photometry of SDSS J1724+5620 with the IAC80 telescope 
at the Observatorio del Teide (Spain) and the 70 cm tele- 
scope of the Astrophysical Institute of Postdam at Babels- 

1 http:/ /www. sdss.org/dr6/dm/flatFiles/spCFrame. html 



3 ORBITAL PERIODS 
3.1 Radial Velocities 

In iRebassa-Mansergas et all l|2007h we measured the ra- 
dial velocities fitting a second order polynomial plus a 
double- Gaussian line profile of fixed separation to the 
Nai AA 8183.27,8194.81 absorption doublet. Free parameters 
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Figure 1. SDSS spectra of the 11 SDSS white dwarf main sequence binaries studied in this work. 



were the amplitude and the width of each Gaussian and the 
radial velocity of the combined doublet. Here we adopt a 
slightly modified approach, using just a single width param- 
eter for both line components. This reduction in the number 
of free parameters increases the robustness of the fits. Radial 
velocities measured in this way for nine WDMS binaries are 
given in Table[2] In addition, we measured the companion 
star radial velocities for SDSS J1724+5620 by means of a 
Gaussian fit to the Ha emission line clearly visible in 22 of 
the SDSS sub-exposures (Sect. [2~Tj) , which are also reported 
in Tabled 

IScargld (|l982l ) periodo grams calculated from the radial 
velocities of each system to investigate the periodic nature 
of the velocity variations contain a number of aliases due to 
the sampling pattern of the observations (Fig. [5] left panel). 
We carried out sine-fits of the form 



Vr = K S , 



2n(t - T ) 



+ 7 



(1) 



to the velocity data sets of each system, where 7 is the sys- 
temic velocity, K Bec is the radial velocity semi-amplitude 
of the companion star, To is the time of inferior conjunc- 
tion of the secondary star, and P or b is the orbital period. 
Several fits were done, adopting the frequencies correspond- 
ing to the strongest peaks in the power spectra as initial 
conditions. The parameters resulting from these fits are re- 
ported in Table ES For SDSS J0052-0053, SDSS J0246+0041, 
SDSSJ0314-0111, SDSS1151-0007, SDSS J1529+0020 and 
SDSS J2339-0020 the sine-fits allow a unique choice of the 
orbital period. For SDSS J1724+5620, the sampling of the 
SDSS spectra is very sparse, resulting in finely structured 
aliases superimposed on a sequence of large aliases spaced 
by Id -1 . A sine fit to the radial velocities started off at 
the 3d -1 alias provides a spectroscopic orbital period of 
7.99243(16) h, which is consistent with the more accurate 
value determined from the photometry (Sect. I3.2|l . The ra- 
dial velocity data for all seven systems folded over their or- 
bital periods are shown in Fig. [2] (right panel). 
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Table 2. Radial velocities measured from the Nal AA 8183.27,8194.81 doublet, except for SDSS J1724+5620, where radial velocities 
measured from the Ha emission line in the 15 min SDSS-subspectra are given. Note that the provided H JD is the resulted from H JD- 
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The periodograms calculated from the radial velocity 
measurements of the PCEB candidates SDSS J0309-0101, 
SDSS J1138-0011, and SDSS J2241+0027 do not reveal any 
significant peak. The low amplitude of the radial velocity 
variations observed in these three objects suggests that they 
may be wide WDMS binaries rather than PCEBs. This issue 
will be discussed further in Sect. 15.31 



3.2 Light curves 

The light curves of SDSS J0314-0111 and SDSS J1724+5620 
display variability with an amplitude of ~ 0.05 mag and ~ 
0.1 mag, respectively. We calculated periodograms for both 
systems using the ORT/TSA command in MIDAS, which 
folds and phase-bins the data using a grid of trial periods 
and fits a series of Fourier t erms to the folded light curve 
l|Schwarzenberg-Czerny|[T996t ) . 

A strong peak is found in the periodogram of 
SDSS J0314-0111 at 3.8 d~\ i.e. 6.3 h (Fig.G3 top left panel). 
Folding the photometric data over that period results in a 
double-humped modulation which we identify as ellipsoidal 
modulation (Fig(3] middle left panel). The detection of ellip- 
soidal modulation indicates that the companion star must 



be filling a significant amount of its Roche-lobe radius and 
have a moderately high inclination, both of which hypoth- 
esis are confirmed below in Sect. [4] The two minima differ 
in depth, as expected for ellipsoidal modulation because of 
the stronger gravity darkening on the hemisphere facing the 
white dwarf. The two maxima are also unequal, which is ob- 
served relatively often in PCEBs, a nd thought to be related 
to the presence of star spots (e.g. iKawka fc Vermes! 120031 ; 
iTappert et~ai1l2007h . We conclude that the strongest peri- 
odicity detected in the photometry of SDSS J0314-0111 is 
consistent with the spectroscopic orbital period (Table|3} ■ As 
the radial velocities span a longer temporal baseline than the 
photometry, they provide the more accurate period measure- 
ment, and we adopt P orb = 6.319 ± 0.015 h for SDSS J0314- 
0111. 

For the analysis of SDSS J1724+5620 we separately 
merged the data for the 7? and I band, normalised both 
data sets to their respective mean values, and finally com- 
bined all data into a single data set. The ORT periodogram 
calculated from the full light curve displays a strong signal 
at 3d -1 , with weak aliases related to the sampling pattern 
inherent to the data (Fig.(3] top right panel). Given the large 
quasi-sinusoidal shape of the modulation, and the high effec- 



6 A. Rebassa-Mansergas et al. 




5 10 

Frequency [d~ 



50 


-50 
200 
100 


-100 

200 
100 


-100 
-200 
200 



-200 

200 
100 


-100 

-200 

100 


-100 

100 


-100 
-200 



?-L,, SDSS J0052-0053 k'H~L ,, ' Porb = 2.74h ' L \ 

- p» t K - 

-. Km ' IS , / 




; 


SDSS J0246+004 


1 t-H ' 1 1 1 ' 1 1 1 tt 

t*-*^* Porb=17.43h 

f ~% s - 




<■ *»•< / 

■ ■■■■■ , ;1 : . ." 


~ 1 1 1 \ 1 

SDSS J0314-01 1 


^f~^-^ Porb=6.32h k 

> S /* 


- • . .. • 


X 


£— N SDSS Jl 151 -000 


r , - ■» * Porb-3.40 h ^ "J 

* - 


: \ r * 




?«■. SDSS J1529 + 0020 ^ -}""■««, Porb = 3.96h sf: 

s ' v .'J 
» ' ■> y - 


L 




**>.. SDSS J1724 + 5620 >**'***- K> Porb = 7.99h V« 

: N >f* x >f* : 

■ x ?. \ ' ■ 


'-- x« 




SDSS J2339-0020 >^">^ Porb=15.74h 

- \ < *s *i — 
\ ** S ** 


\ T 

[ 





0.5 



1 1.5 
Orbital Phase 



Figure 2. Left: Scargle periodograms calculated from the radial velocity variations measured from the Nal absorption doublet in 
SDSS J0052-0053, SDSS J0246+0041, SDSS J0314-0111, SDSS J1151-0007, SDSS J1529+0020, and SDSS J2339-0020, and from the Ha 
emission line in SDSS J1724+5620. The aliases with the highest power are indicated by tick marks. Right: the radial velocity measurements 
folded over the orbital periods of the systems, as determined from the best sine fits (Tablc[3}. 



tive temperature of the white dwarf (Table [3J), we identify 
the observed photometric variability as being due to heat- 
ing effect (or reflection effect as it is often referred to in the 
literature). A sine fit to the combined data (dashed line) 
provided the following photometric ephemeris 



HJD(0 = 0) = 2453865.29(1) + E x 0.3330193(13), 



(2) 



where (f> — refers to the occurrence of minimum light. The 
photometric period of 8h is consistent with the alias pat- 
tern in the periodogram calculated from the Ha radial ve- 
locity variations (Fig. [2] left panel). The phase- folded light 
curve using the photometric period is shown on the middle 
right hand panel of Fig. [3] Based on the above ephemeris, 
the accumulated phase error for the first SDSS subspec- 
trum (Table [2} is 0.024 cycles, which is sufficiently good to 
phase the radial velocities obtained from the SDSS spectra 
with Eq. ((2| . The phase-folded light curve and radial velocity 
curve are shown in Fig. [4] The phase offset between the two 
curves is 0.226±0.005, which is consistent with ~ 0.25 within 
the error of the fit and the ephemeris. A 0.25 phase offset is 
what is expected from the assumption that the photometric 
modulation is related to a heating effect, i.e. maximum light 
and red-to-blue crossing of the radial velocity corresponds 
to the superior conjunction of the secondary star, whereas 
minimum light and blue-to-red crossing of the radial velocity 
corresponds to its inferior conjunction. 

Finally, we show in the bottom panel of Fig. [3] a 3.5 



hours light curve of SDSS J0820+4314. We monitored the 
system through a total of 8.2 hours on three different nights 
(see Table[T]), and found no significant photometric modu- 
lation. This implies that the companion star is significantly 
under-filling its Roche lobe and/or that the binary inclina- 
tion is very low. The orbital period of the system will hence 
need to be measured from a radial velocity study. 



4 BINARY PARAMETERS 



In lRebassa-Mansergas et al.l ((2007) we developed a spectral 
decomposing/fitting technique for the analysis of WDMS bi- 
naries with SDSS spectroscopy. In brief, this analysis deter- 
mines the white dwarf effective temperature, surface gravity, 
mass and radius, the spectral type and radius of the main 
sequence companion, as well as two independent distance 
estimates based on the properties of the two stellar compo- 
nents. Given that the SD SS spectra reduction pipeline ha s 
been improved with DR6 (|Adelman-McCarthv et al.| [2008), 
we have re-analysed here the seven PCEBs for which we were 
able to determine orbital periods. Comparison with t he re - 
sults from DR5 reported in lRebassa-Mansergas et alj (l2007h 
shows that the fit parameters differ slightly, but agree in the 
vast majority of cases within the errors. We report the aver- 
age stellar parameters determined from fits to the multiple 
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Figure 3. Top and middle panels: ORT periodograms and light curves folded on strongest photometric periodicities, P rb = 6.32±0.02h 
for SDSS J0314-0111 and P orb = 7.99246 ± 0.00003 h for SDSS J1724+5620. Bottom panel: the light curve of SDSS J0820+4314 is 
essentially flat, and a periodogram calculated from this data does not contain any significant periodicity. 



SDSS spectra in Tabled Rewriting Kepler's third law (as- 
suming common notation), 



(M wd sin z) 
(M wd + M sc . 



P uK* 
2-kG 



iVs, 



M wd \2^GJ 



1/3 



(M wd + Ms, 



n2/3 



(3) 



(4) 



and adopting the orbital period determined from the analy- 
sis of the radial velocities and the photometry, as well as the 
masses from the analysis of the SDSS spectra, we are now 
able to estimate the orbital inclinations for the seven sys- 
tems in Tabled Two systems require some additional notes: 
SDSS J0314-0111 contains a DC white dwarf, and hence no 
white dwarf parameters could be determined from the spec- 
tral analysis, and we assume M w d = O.65M0 for the esti- 
mate of the inclination. In SDSS J1724+5620 the hot white 
dwarf is irradiating the companion star, and the implications 
are discussed in more detail in Sect . 15.21 

Inspecting Eq.Q, it is clear that the dominant uncer- 
tainties in the inclination estimates are only M w d and M sec , 
as the orbital periods and K scc velocities are accurately de- 
termined. The primary uncertainty in the inclination esti- 
mates is M scc , as it is based on the spectral type of the 
companion sta r, adopting the spectral type-ma ss-radius re- 
lation given by |Rebassa-M anserga s~et al.l l|2007h . M w d is de- 
termined from fitting the Balmer line profiles, and relatively 



well constrained. Given that sin i oc Ms/c , but sin i oc M w d, 
the relative weight of the uncertainty in M scc is alleviated. 
We estimate the uncertainties on the binary inclinations in 
Table|4] by assuming in Eq.Q the range in M scc implied 
by a spectral type uncertainty of ±0.5 spectral classes plus 
the associated errors in the spectral type-mass relation, as 
well as the range in M w d resulting from the Balmer line 
profile fits. Given the high estimate for the binary inclina- 
tion, SDSS J1529+0020 is a good candidate for photometric 
follow-up observations probing for eclipses in its light curve. 

In addition, knowing that M ae c/M w d = ^wd/^sec = q, 
we can also estimate the expected orbital velocity of the 
white dwarf K v d (Tabled . The predicted K^d amplitudes 
could easily be measured, e.g. using HST/COS observations, 
and such measurements would be very valuable to improve 
the overall constraints on the system parameters of these 
PCEBs. Finally, we can estimate the orbital separations and 
Roche lobe radii of the se condary stars from Kepler's third 
law and Eggleton's l|l983T ) expression 



(0.49 q 2/A 



0.6 g 2 / 3 + ln(l + 9 1 



/■■n 



(•») 
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Table 4. Average binary parameters obtained for th e seven PCEBs which have orbita l period and K BCC measurements. M w( j, M Bec , R Be c, 
spectral type, T e ff and logg are obtained following Rebassa-Mansergas ct al. (2007), except for SDSS J0314— 0111, where we assume a 
WD mass of 0.65 ± 0.1 Mq (see text for details). P or b and K scc are measured in Sect.[3]of this work. Estimates of orbital separations 
a, q, /f w(i , secondary Roche lobe radius Rl soc and inclinations are obtained from the equations given in Sect. [4] SDSS J1724+5620 is a 
particular case, as the inner hemisphere of the companion is irradiated. Constraints on its orbital parameters are obtained assuming a 
spectral type between M3-5, and considering different K scc corr values for each mass and radius (see Sect. l5.2ll . An additional constraint 
for its inclination comes from the fact that SDSS J1724+5620 is not eclipsing. 
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Table 3. Orbital periods, semi-amplitudes K sec , systemic veloc- 
ities 7 sec , and reduced \ 2 from sine fits to the Na I doublet radial 
velocity data for the strongest two to three aliases in the peri- 
odograms shown in Fig. [2] The best-fit values are set in bold. 
The SDSS Ha radial velocities of SDSS J1724+5620 are folded 
on the photometric orbital period obtained in Sect. 13.21 (which is 
more accurate than the spectroscopically determined value of the 
orbital period). Note also that the semi-amplitud measurement 
of this system is underestimated, as it comes from RV measure- 
ments of Ha emission from the irradiated face of the companion 
(see Sect. [52}. 



5 NOTES ON INDIVIDUAL SYSTEMS 

5.1 SDSS J0052-0053, a detached CV in the 
period gap? 

SDSS J0052-0053 has the shortest orbital period, the small- 
est radial velocity amplitude, and lowest inclination in 
our sample (Fig. [2] Table [IJ. Another intriguing feature of 
SDSS J0052-0053 is that the Roche lobe secondary radius 
and the secondary star radius overlap within the errors (see 
Table|4]). The SDSS and Magellan spectra certainly rule out 
ongoing mass transfer, i.e. that SDSS J0052-0053 is a dis- 
_guised cataclysmic variable (CV). Two possible scenarios 
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1.2 



^eould apply to the system. SDSS J0052-0053 could be ei- 
ther be a pre-CV that is close to develop in a semi-detached 
Configuration, or it could be a detached CV in the 2-3 h pe- 
3-74i oc i g a p. Standard evolution m odels based on the disrupted 

■ 5 fti agnetic braking scenar io (e.g. iRappaport et al.ll 19831 ; iKolbl 

1 ' ^9931 ; iHowell et al.|[200ll ) predict that CVs stop mass trans- 
' ,er and the secondary star shrinks below its Roche lobe ra- 

22 jams once they evolve down to P or b — 3h. Subsequently, 
27 Jhese detached CVs evolve through the period gap, until the 
l.lgompanion star fills its Roche lobe again at P or b — 2 h. Pop- 
34. elation models based on the disrupted magnetic braking hy- 
3.5lj»othesis predict that the ratio of detached CVs to pre-CVs 
64.^ w ith appropriate companion star masses to (re-)initiate 
2 6-niass transfer at P or b — 2 h) should be > 4 (|Davis et al.l 
hence a substantial number of such systems is ex- 
^'?>ected to exist. So far, only one ot her system with simil ar 
°'g properties is known, HS 2237+8154 |Gansicke et alj|2004) . 

0.69 



5.2 



SDSS J1724+5620, a PCEB with a strong 
heating effect 



SDSS J1724+5620 contains the hottest white dwarf among 
our sample of PCEBs, which affects the determination of 
its system parameters in several ways. Firstly, irradiation 
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Figure 4. Top panel: The IAC 80 and AIP 70 cm photome- 
try of SDSS J1724+5620 folded over the photometric period of 
7.9924632 ± 0.0000312 h. Middle pane: the equivalent width vari- 
ation of the Ha emission line measured from the SDSS spectra of 
SDSS J1724+5620 (Table folded over the photometric period. 
Maximum equivalent width occur roughly in phase with the max- 
imum in the light curve (see Sect. 15.21 for details). Bottom panel: 
the radial velocity variation of the Ha emission line. The rela- 
tive phasing with respect to the photometry is consistent with 
the photometric modulation being caused by irradiation of the 
secondary star by the hot white dwarf. 



heating the inner hemisphere of the companion will cause 
it to appear of earlier spectral type than an unheated star 
of same mass. This effect is observed in the analysis of the 
three SDSS spectra available in DR6. The spectral decompo- 
sition types the companion star as M3-4, and combining the 
flux scaling factor of the M star template and the spec tral 
type-radius relation (see iRebassa-Mansergas et al.ll2007j . for 
details) indicates an average distance of d soc = 633 ± 144 pc. 
This is roughly twice the distance implied from the model 
fit to the residual white dwarf spectrum, d w d = 354 ± 15 pc, 
indicating that the spectral type of the companion star de- 
termined from the SDSS spectrum is too early for its actual 
mass and radius, as expected for being heated by the white 
dwarf. Fixing the spectral type of the companion to M5, 
the spectral decomposition results in d acc = 330 ± 150 pc, 
consistent with the distance based on the white dwarf fit. 

An additional complication is that the radial velocity of 
the companion star was measured from the Ha emission line, 
as the Na i absorption doublet was too weak in the individual 
900 s SDSS spectra. The equivalent width of the Ha emission 
shows a noticeable variation as a function of binary phase, 
with maximum equivalent with near the maximum in the 
light curve, indicating that the Ha emission is concentrated 
on the inner hemisphere of the companion star (Fig. [4}. A 
phase offset of 0.046T0.008 is observed between the photom- 
etry and the Ha equivalent width, which we attribute to sys- 
tematic problems in measuring the equivalent width of the 
Ha embedded in the photospheric absorption line from the 



white dwarf, given the poor quality of the individual SDSS 
subspectra. As a consequence of Ha predominantly origi- 
nating on the inner hemisphere of the companion star, the 
observed K BCC is an underestimate of the true radial veloc- 



ity amplitude (e.g. IWade & Home! 1988 


: Orosz et al.lll999l: 


IVennes et al.l Il999l: lAungweroiwit et al. 


120071). The radial 



velocity amplitud e of the secondary st ar's centre of mass 
can be written as l|Wade fc Hornelll988h : 



K, 



scc,corr 



(l + q)(AR/a) 



(6) 



where AR is the displacement of the centre of light from 
the centre of mass. AR — implies that the centre of light 
and the centre of mass coincide, whilst AR = i? soc gives 
the maximum possible displacement. If one assumes that 
the irradiated emission on the secondary is distributed uni- 
formly over the inner hemisphere, and that the contribu- 
tion of the irradiati on is zero on its un-irradiated face, then 
AR = (4/37r)fl sec i|Wade fc Home! 1 19881 ; lOrosz et al.lll999l : 
IVennes et alJll99sh . Assuming that the spectral type is in 
the range M3-5, different combinations of secondary mass 
and radius, and radial velocity amplitude, K SCCiCorr , can then 
constrain the orbital parameters of SDSS J1724+5620 (Ta- 
ble©. 



5.3 SDSS J0309-0101, SDSS J1138-0011 and 

SDSS J2241+0027: wide WDMS binaries? 

SDSSJ0309-0101, SDSSJ1138-0011 and 

SP SS J2241+0027 were fla g ged as PCEB candidates 
by IRebassa-Mansergas et ail (|2007t ) on the base of a 3a 
radial velocity variation in between their different SDSS 
spectra, as measured from either the Ha emission line or the 
Nai absorption doublet from the companion star. However, 
the additional intermediate-resolution spectra taken for 
these three objects (Table [TJ do not show a significant 
radial velocity variation (Table[2] Fig. [5}. It is therefore 
important to review the criterion that we used in our first 
paper to identify PCEBs from repeated SDSS spectroscopy 
in the light of two subtle changes. 

On the one hand, we have slightly modified the pro- 
cedure to fit the Nai absorption doublet, as outlined in 
Sect. 13. 1| By comparing t he Na i DR5 radial velocities in 
IRebassa-Mansergas et al] (|2007l ) with those obtained with 
the new procedure for the same spectra, we find an aver- 
age difference of Skms" 1 with a maximum difference of 
10.5 kms^ 1 . In all cases the measurements agree within the 
errors. 

On the other hand, we w ere u sing DR5 spec- 
tra in IRebassa-Mansergas et all l|2007l ). but the anal- 
ysis carried out here was done using DR6 spectra, 
which were processed with a d ifferent reduction pipeline 
|Adelman-McCarthv et~aH 120081 '). A comparison between 
the DR5 and DR6 rad i al ve locity values obtained in 
IRebassa-Mansergas et ah! (|2007l ) and in this work respec- 
tively (both measured following our new procedure) provides 
in this case an average difference of 6.5 km s -1 , with a maxi- 
mum of 22kms-\ A gain, the radial velocity measurements 
agree, with the exception of a single spectrum, within the 
errors. 

The conclusion from comparing our two methods, and 
the two SDSS data releases, is that the radial velocity 
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Table 5. Parameters of the PCEBs derived from calculating their post-CE evolution according to Schrcibcr & Gansickc (2003). t CO ol and 
t s d are the cooling age and the predicted time until the system enters the semi-detached CV phase respectively. P st j denoted the zero-age 
CV orbital period while Pcb is the orbital period the system had when it left the CE phase. As discussed in the text and displayed 
in Fi g[6l the predicted evolution depends sensitively on the mass and the radius of the secondary star, in particular as the errors of all 
estimated secondary masses overlap with the value assumed for the fully convective mass limit. We therefore give values for the both 
scenarios. 











CMB 




GR 


Name 


-^"orb 


^cool 




Pce 




Psd 






SDSSJ 


[days] 


[years] 


[days] [days] 


[years] 


[days] [days] 


[years] 


0052-0053 


0.114 


4.2 x10 s 


0.125 


0.46 


0.0 


0.114 


0.145 


0.0 


0246+0041 


0.726 


3.1x10 s 


0.146 


0.78 


l.lxlO 9 


0.113 


0.73 


8.1xl0 10 


0314-0111 


0.263 




0.12 




3.7xl0 7 


0.11 




6.2xl0 9 


1151-0007 


0.142 


5.0x10 s 


0.12 


0.55 


2.7xl0 6 


0.07 


0.15 


1.8xl0 9 


1529+0020 


0.165 


1.3x10 s 


0.11 


0.41 


4.7xl0 6 


0.10 


0.17 


2.5xl0 9 


1724+5620 


0.333 


3.2xl0 6 


0.12 


0.34 


6.1xl0 7 


0.11 


0.33 


1.7xl0 10 


2339-0020 


0.656 


5.1x10 s 


0.12 


0.74 


1.0x10 s 


0.11 


0.66 


6.5xl0 10 



measurements obtained are in general consistent within 
their errors. However, the offsets can be sufficient to move 
a given system either way across our criterion to iden- 
tify PCEB candidates, being defined a 3a radial varia- 
tion between their SDSS spectra. This is specifically the 
case for systems with either low-amplitude radial velocity 
variations, or faint systems with noisy spectra. An addi- 
tional note concerns the use of the Ha line as probe for 
radial velocity variations. iRebassa-Mansergas et ail (|2007T ) 
identified SDSS J2241+0027 as a PCEB candidate on the 
basis of a large change in the Ha radial velocity be- 
tween the two available SDSS spectra. However, the veloc- 
ities obtained from the Nai doublet did not differ signifi- 
cantly. Inspecting the spectra a gain confirms the results of 
IRebassa-Mansergas et all |2007h . This suggests that the ab- 
sorption lines from the secondary star are a more robust 
probe of its radial velocity. 

A somewhat speculative explanation for the shifts found 
from Ha radial velocity measurements is that the main se- 
quence star is relatively rapidly rotating, and that the Ha 
emission is patchy over its surface. A nice example of a 
this effect is the rapidly rotating (P — 0.459 d) active M- 
dwarf EY Dra, which displays Ha radial velocit y variation s 
with a peak-to-peak amplitude of ~ 100 km s -1 (|Eibell 19981 ). 
In order to explain the Ha radial velocity shifts of a few 
lOkms" 1 observed in e.g. SDSS J2241+0027, the compan- 
io n star should be rotating with a period of ~ 1 d. According 
to ICardini"fc Cassatellal (2007), low-mass stars such as the 
companions in our PCEBs that are a few Gyr old are ex- 
pected to have rotational periods of the order of tens of days, 
which is far too long to cause a significant Ha radial veloc- 
ity variat ion at the spectral res olution of the SDSS spectra. 
However, ICardini fc Cassatellal 's study was based on single 
stars, and the rotation rates of the companion stars in (wide) 
WDMS binaries is not well established. 



6 PCEB EVOLUTION 

Following ISchreiber fc Gansickd l|2003h we determine the 
cooling age and the future evolution of the new PCEBs. 
We assume classical and disrupted magnetic braking ac- 
cording to the standard theory of CV evolution (e.g. 



SDSSJ0309-0101 SDSSJ1 138-0011 SDSSJ2241 +0027 

100 : ! 



50 - 1 -' 

■ 100 o o 

Figure 5. SDSS radial velocities (triangles) along with the Nai 
radial velocities measured in this work (solid dots) of SDSS J0309- 
0101, SDSS J1138-0011 and SDSS J2241+0027. The data at hand 
suggest that this three systems are wide WDMS binaries instead 
of PCEBs. 

Verbunt fc Zwaanl Il98ll ). We used the cooling tracks of 
Wood! l| 19951 ) and find that most of the PCEBs left the 
CE about 1 — 5 x 10 s years ago, the only exception be- 
ing SDSS J1724+5620 which appears to be much younger. 
The time still needed to enter the semi-detached CV config- 
uration is shorter than the Hubble time for all the systems, 
qualifying them as representatives for the progenitors of the 
current CV population. The numbers we obtained from the 
theoretical analysis are summarised in Table[5] The esti- 
mated secondary masses are generally very close to the fully 
convective boundary (for whic h we assume M cc = 0.3Ma, 
see however the discussi on by iMullan fc MacDonaldl l200ll 
and lChabrier et al. I l2007l on how this mass limit may be af- 
fected by magnetic activity) and the uncertainties involved 
in the secondary mass determination are quite large. The 
predicted evolution of all systems should therefore be con- 
sidered highly uncertain as it is not clear whether magnetic 
braking applies or not. We therefore give in Table[5]the val- 
ues obtained for both scenarios. The effect of not know- 
ing whether magnetic braking or only gravitational radia- 
tion will drive the evolution of the systems is illustrated 
in Fig. HI] which shows the expected post CE evolution for 
SDSS J1152-0007 and SDSS J0246+0041 for both cases. 

Inspecting Tabled the case of SDSS J0052-0053 is par- 
ticularly interesting. The estimated secondary radius is con- 
sistent with its Roche-lobe radius and the system is sup- 
posed to be very close to the onset of mass transfer, which 
is reflected by t s( j = within the accuracy of our calcula- 
tions. Clearly, SDSS J0052-0053 may be a detached CV in 
the period gap, or a pre-CV that has almost completed its 
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PCEB lifetime. To estimate the probability for SDSS J0052- 
0053 being either a period-gap CV or a PCEB close to en- 
tering for the first time a semi-detached state, we assume 
a steady state binary population. In that case, the number 
of detached CVs in the gap sh ould be rou ghly equal to the 
number of CVs above the gap (|Kolblll993r i. In addition, the 
number of PCEBs in the gap should be roughly equal to the 
number of accr eting CVs i n the gap. According to CV pop- 
ulation studies (Kolb 1993), the number of detached CVs in 
the gap is approximately five times higher than the number 
of CVs in the gap and the probability for SDSS J0052-0053 
being a detached CV in the gap rather than a PCEB is 
~ 80%. This result is a simple consequence of the fact that 
all CVs with donor stars earlier than ~M3 will become de- 
tached CVs in the gap while only a rather small fraction 
of PCEBs, i.e. those with secondary stars of spectral type 
~M3.5-M4.5, produces CVs starting mass transfer in the 
period gap. This is obviou sly only a rough e stimate. A de- 
tailed population study bv lDavis et al.l (|2008f ) confirms that 
the ratio of detached CVs to pre-CVs within the CV pe- 
riod gap is ~ 4 — 13, depending on different assumptions of 
the common envelope ejection efficiency, initial mass ratio 
distributions, and magnetic braking laws. 



7 DISCUSSION 

As outlined in the introduction, despite some recent progress 
our understanding of th e CE phase is st ill very limited. The 
classical a prescription (|Webbink| |l984|) that is based on a 
simple energy equation seems to be unable to explain the 
observed existence of double white dwarfs with stellar com- 
ponents of_com^arable masses and rather long orbital peri- 
ods fsee lWebbin k 2007, for more details). This disagreement 
motivated iNelemans et al.l (|200pJ) to develop an alternative 
prescription based on the conservation of angular momen- 
tum, the so-called 7 algorithm. Binary population synthesis 
codes based on the a prescription as well as those assum- 
ing the 7 algorithm predict the existence of a significant 
number of PCEBs with orbital periods longer than a day. 
Using the a mechanism, IWillems fc Kolbl (|2004| ) calculated 
the expected period distribution of the present-day WDMS 
binaries in the Galaxy at the start of the WDMS binary 
phase. Their Fig. 10 clearly shows that the predicted PCEB 
distribution peaks around P or b ~ 1 day, but also has a long 
tail of systems with up to ~ 100 d. The 7 algorithm, on the 
other hand, predicts an increase of the number of PCEBs 
with increasing orbital period up to P rb~ 100 days, i.e. the 
7 prescription predicts the existence of e ven more long or- 
bital period (P orb > 1 day) PCEBs (see INelemans fc Tout! 
120051 ; (Maxted et aI1l2007h . 

We have measured a total number of nine orbital pe- 
rio ds along this work a nd the recently published paper 
by ISchreiber et all l|2008h . All of them have short orbital 
periods of less than a day. Radial velocity variations are 
much easier to identify in short orbital period systems, so 
our sample is expected to be biased towards shorter or- 
bital periods. The crucial question is whether our obser- 
vational finding, the paucity of PCEBs with P or b > 1 d, is 
whether the result of that bias, or whether it reflects an 
intrinsic feature of the PCEB population. To answer this 
question we performed the following analysis. First, we car- 




Figure 6. The predicted evolution of two PCEBs. The orbital 
period gap observed in the orbital period distribution of CVs 
is shaded. As the secondary masses of SDSS J1151-0007 and 
SDSS J0246+0041 are close to the fully convective boundary, 
we calculated the evolution assuming classical magnetic braking 
(CMB) and assuming only gravitational radiation (GR). Obvi- 
ously, the calculated orbital periods at the end of the CE phase 
(Pce), the evolutionary time scale, and the expected zero-age CV 
orbital period (-P s d) differ significantly. 



rie d out Monte-Ca rl o sim ulations similar to those presented 
in ISchreiber et al.l |2008) but assuming a resolution of 
15km/s, corresponding to the typical error in the radial ve- 
locity measurements from SDSS spectra. Clearly, as shown 
in Fig. [7] the bias towards short orbital period systems is 
larger due to the lower resolution of the SDSS spectra and, 
as a consequence, radial velocity variations in multiple SDSS 
spectra of systems with P or b > 10 days are more difficult 
to detect. However, the probability to detect 3cr radial ve- 
locity variations for systems in the orbital period range of 
Porb ~ 1 — 10 days still is ~ 20 — 40%. Second, we assumed 
a uniform orbital period distribution and integrated the 3a- 
detection probability (lowest dotted line) for systems with 
Porb ^ 1 day and those with 10 days > P or b > 1 day. In 
other words, we assume that there is no decrease for 10 days 
Porb > 1 day and calculate the total detection probability 
for 10 days > P rb > 1 day and P or b ^ 1 day. We find that 
if there was indeed no decrease, the fraction of PCEBs with 



2 In [Schrcibcr et al.l (|200d . their Fig. 8) we underestimated the 
detection probability. Therefore, the correct detection probability 
for a resolution of 6km/s as given in Fig. [7] is higher than our 
previous estimates. 
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Figure 7. Monte-Carlo simulations of the detection probability 
of significant radial velocity variations assuming measurement ac- 
curacy of 6 kms -1 , corresponding to our previous VLT/FORS ob- 
servations llSchreib er et al .ll2008t) 15km/s, as a ppropriate for the 
SDSS spectra llRebass'a-Mansereas et alj|2007h . The three lines 
correspond to 1, 2, or 3<r significance of the radial velocity vari- 
ation. Clearly, even the PCEB identification based on multiple 
SDSS spectra should be sensitive to PCEBs with orbital periods 
of ~ 1 — 10 days. 



Parb > 1 day in our sample should be ~ 84% and the number 
of PCEBs with Porb > 1 day among our nine systems should 
be ~ 7.6 ±2.8. The result of our observations, i.e. no system 
with P or b > 1 day among nine PCEBs disagrees with the 
hypothesis (i.e. there is no decrease) by 2.7a. This indicates 
that the measured lack of PCEBs with P or b > 1 day might 
indeed be a feature of the intrinsic population of PCEBs. 

How do our results for PCEBs from SDSS compare with 
the sample of all known PCEBs? Figure|S] shows in gray the 
orbital period distribution of all 41 known PCEBs consist- 
ing of a M or K-star plus a white dwarf from the latest 
version of the iRitter fc Kolbl i|2003T ) catalogue (V7.9), in- 
cluding also the seven new periods from this pap er, and 
the new recent discovery by ISteinfadt et al.l |2008h . As we 
are interested in analysing the bias of the sample towards 
shorter orbital periods in the context of detection proba- 
bility, we superimpose in black a subsample of ten SDSS 
PCEBs, identified from radial velocity snapshots of SDSS 
WDMS, representing a homogenous selec tion mechanism 
(seve n from the current paper, two from ISchreiber et al] 
20081 . one [SDSS J112909. 50+663704.4] from lRavmond et all 
2003). So far, the period distributions of the two subsets 
seem very similar, with a steep decrease of the number of 
systems at P or b ~ 1 day. The total sample of PCEBs con- 
tains only six systems with orbital periods larger than one 
day (~15% of the entire sample). This suggests that the 
number of PCEBs decreases for P or b > 1 day, implying that 
the 7 algorithm in its present form might not be an adequate 
description of the CE phase, and that the common envelope 
efficiency in the energy equation used in the a prescription is 
perhaps smaller than assumed previously (lower efficiencies 
lead to stronger shrinkage of binary separations). However, 
one should bear in mind that the previously known PCEBs 
have been discovered by various met hods, and therefore do 
not form a representative sample (see lSchreiber fc Gansickel 




log(Porb) [d] 



Figur e 8. The pe riod distribution of 41 PCEBs from the 
IRitter &: Ko lb (2003) (V7.9) catalogue containing a white dwarf 
plus M or K companion star (gray). The 7 systems studied in this 
work and the recent P CEB discovery SDSS J143547.87+373338.5 
iSteinfadt et alj|2008j ) are also included. In black is shown a sub- 
sample of SDSS PCEBs obtained from RV variation studies. Both 
subsamplcs illustrate a clear lack of systems at orbital periods 
greater than 1 day. The 7 new systems fr om this work and the 
the two systems in ISchreiber et al.l |2008() are indicated as tick 
marks in the top of the figure. 

l2003t for details). The sample of PCEBs selected in a homo- 
geneous way from the SDSS is still small, and the current 
sample of WDM S binaries from SDSS involves some selec- 
tion effects, too (|Schreiber et alj|2007f ). impeding a definite 
conclusion at this point. However, work is underway to en- 
large the parameter space of the SDSS W DMS binaries in 
term s of ages and secondary star masses (Schrcibcr et al. 
2007), to increase the number of PCEBs with orbital period 
measurements, and to model the remaining selection effects, 
paving the way for a more quantitative assessment of CE 
and post-CE evolution. 



8 CONCLUSIONS 

We have presented a study of 11 PCEB previously identified 
candidates, and measured the orbital periods of six and one 
of them from their Na i doublet radial velocity variations and 
its differential photometry, respectively. Combining the K scc 
velocity amplitudes with the results from spectroscopic de- 
composition/fitting of their SDSS spectra, we constrained 
the binary parameters of the seven systems for which we 
could determine orbital periods. No radial velocity varia- 
tions were detected for three PCEB candidates, suggesting 
that they may be w ide binaries. We revisited the P CEB 
candidate selection of lRebassa-Mansergas et al.l (|2007t ). and 
conclude that candidates with low amplitude velocity vari- 
ations, noisy SDSS spectra, and radial velocity shifts that 
only show up in the Ha emission line definitely need addi- 
tional follow-up spectroscopy for the confirmation/rejection 
of their PCEB nature. Finally, we have had a first look at the 
period distribution of PCEBs from SDSS identified from ra- 
dial velocity snapshots of SDSS WDMS binaries, and noted 
that none of the ten systems published so far have an orbital 
period > 1 d. Using a Monte-Carlo simulation, we demon- 
strated that our method of finding PCEBs should be effi- 
cient also for systems with periods > 1 d, and that, sub- 
ject to small number statistics, it appears that the PCEB 
period distribution peaks at P or b < 1 d. This is in agree- 
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ment with the period distribution of the previously known 
PCEBs, which represents a heterogenous mix of systems 
identified by various different methods. In contrast, current 
binary population models predict a large number of PCEBs 
with orbital periods > 1 day in clear disagreement with the 
currently observed sample. Additional effort is needed to im- 
prove the size of the sample of known PCEBs, as well as to 
fully model its selection effects, but with more than fOOO 
SDSS WDMS binaries to draw from, the outlook for more 
quantitative tests of CE evolution seem promising. 
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